Dislocations and vortices in pair density wave superconductors 
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With the ground breaking work of the Fulde, Ferell, Larkin, and Ovchinnikov (FFLO), it was 
realized that superconducting order can also break translational invariance; leading to a phase 
in which the Cooper pairs develop a coherent periodic spatially oscillating structure. Such pair 
density wave (PDW) superconductivity has become relevant in a diverse range of systems, including 
cuprates, organic superconductors, heavy fermion superconductors, cold atoms, and high density 
■ quark matter. Here we show that, in addition to charge density wave (CDW) order, there are PDW 

ground states that induce spin density wave (SDW) order when there is no applied magnetic field. 
Furthermore, we show that PDW phases support topological defects that combine dislocations in 
£S| ■ the induced CDW/SDW order with a fractional vortex in the usual superconducting order. These 

defects provide a mechanism for fluctuation driven non-superconducting CDW/SDW phases and 
conventional vortices with CDW/SDW order in the core. 



PDW superconductivity, of the kind originally discussed by FFLO 0, Q , is believed to exist in the heavy fermion 
superconductor CeCoIns uA Nfl and in the organic superconductor k-(BEDT-TTF) 2 Cu(NCS) 2 0,0. It is also believed 
be relevant in cold atoms [3, Hi and in the formation of color superconductivity in high density quark matter Also, 
^ . PDW superconductivity with quite a different origin has been found in microscopic theories of correlated electronic 
materials 13, 11. 12 1. Such order is believed by some to be competing with conventional <i-wave superconductivity in 
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. underdoped high temperature cuprate superconductors. Specifically, PDW order may appear as an alternate ground 
state near a hole doping of 1/8, where some sort of charge order is often observed lH, 1J]. Given its relevance, it is 



important and useful to address the properties of PDW order from a phcnomcnological point of view. That is the 
primary goal of this work. In particular, we address two questions that apply to both FFLO superconductors and 
'-^J . the underdoped cuprates: 1) What symmetries are broken by PDW order? and 2) What are the properties of the 
vortex-like topological defects in the PDW ordered phases? The answer to the first question reveals that CDW or 
SDW order often must accompany the PDW order and provide a means to identify the PDW phase. The answer 
to the second question turns out to be non-trivial. It is clear that there will be superconducting vortices and these 
have been studied in the past (see Ref. [6\ for an overview). However, it is also well known that the periodic order 
has dislocations as natural topological defects. This has not been given much consideration in the context of PDW 
order. Here we show that the topological defects of PDW superconductors contain not just vortices and dislocations, 
but also combinations of fractional vortices and fractional dislocations. We argue that these fractional defects play an 
important role when considering fluctuations in two-dimensions (for example, leading to non-superconducting CDW or 
SDW phases) and that they play an important role when considering the physics of the usual superconducting vortices 
^vq (for example, leading to the appearance of CDW or SDW order inside the vortex core). Finally, we address an issue 
that is mainly of relevance to the underdoped cuprates. In particular, we examine the role of competition between 
0^ . PDW order and translational invariant d-wave superconductivity. We show that this competition preferentially selects 
PDW phases with CDW order. Furthermore, in addition to the CDW order stemming from the PDW order parameter, 
^ . the coexistence of PDW and d-wave superconductivity leads to either SDW or additional CDW order. 

While our main results apply more generally to PDW superconductors, for concreteness we consider an example 
motivated by theoretical and experimental proposals for the underdoped cuprates that exhibit some type of charge 
order. In particular, we consider a PDW superconductor in three dimensional tetragonal system with a lattice spacing 
a for the two-dimensional square lattices. The PDW order is taken to be either commensurate, with a periodicity of 
Tol 11 , Tj^, 14 1 , or incommensurate (this also allows the theory to describe a variety of FFLO phases found to be 



stable in two-dimensions |6|). Furthermore, we will take this order to be aligned along the x (taken to be along a 
two-fold symmetry axis) or equivalent directions. The PDW order parameter is written as {Aq x , Aq v , A_q x , A_Q y } 
describing PDW order with corresponding wavevectors {Q x , Q y , — Q x , — Q y }- The Ginzburg Landau free energy 
density is constructed by imposing translational, gauge, time-reversal, parity, and tetragonal rotation symmetries. 
This yields 

/= «EJA Q J 2 + /3 1 (E l A Q J 2 ) 2 + /3 2 E l<J |A Q J 2 |A Q J 2 + /3 3 (|A Q J 2 |A_ Q J 2 + |A Q J 2 |A_ Q J 2 ) (1) 
+/3 4 [A Qx A_ Qx (A Q! A_ Qy )* + (Aq x A_qJ* A Qy A_qJ 
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where (— I) 1 is —1 for = : Q x and 1 for = ±Q„, D = — zV — 2eA, and B = V x A. The difference between the 
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TABLE I: Properties of PDW Ground States. All possible PDW ground states and accompanying CDW and SDW order. 
The third column shows the parameter regions for which these phases are stable. In the fourth and fifth columns other modes 
can be found by using the relationships pq = (p_q)* and Sq = (SIq)* . The fifth column shows the minimum flux contained 
by a topological defect. In Phases 4 and 5, the $o/4 defects experience a confinement potential but can exist at short length 
scales and may consequently become physically relevant. 



commensurate and incommensurate cases arises at eighth order. This implies that sufficiently near the mean field T Cl 
this difference can be ignored. However, there are situations where this difference can be important and this will be 
discussed later (in particular, when considering vortex-type excitations in a two-dimensional PDW superconductor). 

An important property of the above free energy is that it contains a U(l) x U(l) x U(l) symmetry. It is this feature 
that gives rise to vortices and dislocations as well as the appearance of the combined fractional vortex and fractional 
dislocations mentioned above. Below, we will provide a physical picture for the topological defects of Eq. Q] (and also 



for the related half-flux vortices found in Ref. 15|). Prior to understanding these defects, it is important to first 
understand the possible PDW ground states and also useful to understand the induced charge density wave (CDW) 
and spin density wave (SDW) order that arises. To complete the latter goal, we include the free energy for CDW and 
SDW order and their coupling to the PDW order. These are also determined by the same symmetry requirements as 
above. The contribution to the free energy from the CDW order, pq is 

E^ a Q,<°Q*<°-Q' +< J / 9 Q 3 -Q I [ A Q I ( A Q 3 )* + a -qA A -Q,T]- ( 2 ) 

For our choice of PDW order, only SDW order with moments oriented along z can be induced. Denoting the SDW 
order as Sq, the relevant free energy is 

^{^•a s Qi 54^ Qi +< i 54_ Qi i[A Qi (A Qi r- A_ Qj (A_ Qi r]}. (3) 

Throughout this work we consider the case that the CDW and SDW order is induced by the PDW order (qq , ctq > 0) 
and thus is representative of the symmetries broken by the PDW phase transition. In this case they are given by 

a p - 

PQ*-Q* = ~2^— [^WW + A_ Qj .(A_ Q< )*] (4) 



= -^F^[Aq,(Aq,)* - A_ Q ,(A_ Q .)*]. (5) 

Eqs. H] and [5] reveal that the phase difference between two different components Aq can be interpreted as the phase 
of either the CDW or SDW order. In general, it is possible that there is also an intrinsic CDW, or SDW order (for 
which a P Q < or «q < is possible). This is also of interest but is not considered in this work. 

Remarkably, it is possible to find all the possible PDW ground states of the general free energy defined by Eqs. Q] 
(this is not true, for example, for a general phenomenological theory of superfluid 3 He [3]). These ground states and 
the induced CDW and SDW orders are listed in Table |TJ Note that we have listed the lowest Fourier components 
of the induced SDW and CDW order, higher Fourier harmonics will in general appear but these should be smaller 



3 



in magnitude. In principle, the induced CDW and SDW order can be measured and therefore used to identify what 
type of PDW order appears. One result that has not previously been highlighted is the existence of PDW phases 
that break time-reversal symmetry. In particular, Phases 2 and 5 exhibit this behavior and the broken time reversal 
symmetry is readily apparent through the associated SDW order. At the mean field level, this theory predicts that 
superconductivity appears at the same transition temperature as the SDW and CDW order. However, fluctuations 
can separate these transitions and a specific mechanism for this will be discussed later. It is interesting that related 
SDW and CDW order has been observed in a variety of cuprates 17 , 18, 19[ and a systematic experimental study of 



such order in a single material would prove useful to identify any PDW order parameter. 

We now address the nature of the competition between translational invariant d-wave superconductivity and PDW 
superconductivity relevant only in the context of the cuprates. The lowest order coupling between these two order 
parameters is given by the following free energy density: 

/ c = /3 cl ^|A d | 2 |A Q J 2 + / 3 c2 [A/(A Q ,A_ Q . E +A Qy A_ Q J* + (A/)*(A Qi A_ Qi .+A QB A_ Qy )]. (6) 

i 

The relevant feature of this coupling is that the term with coefficient /3 C 2 can always be made negative by an appropriate 
choice of phases of the different order parameters. Consequently, any PDW phase for which this coupling is non-zero 
will lower its energy through this coupling. Only Phases 3 and 4 above lead to a non-zero (32c coupling, these are the 
two phases with PDW induced CDW order. When PDW and <i-wave superconductivity co-exist, then the interplay 
between these two order parameters will lead to CDW or SDW order that appears in addition to the CDW order that 
already exists in Phases 3 or 4. In particular, if the coefficient (3 C 2 > then the additional order will be SDW order at 
the wavevectors of the PDW order; if (3 C 2 < 0, then there will be no induced SDW order and additional CDW order 
will appear at the wavevectors of the PDW order. 

With the understanding of the possible PDW mean-field ground states, we now turn to understanding the topological 
defects. Single valuedness of the wavefunctions implies that these can be found by allowing the phases </>i, </>2j or </>3 
in Table U to have an integral multiple of phase winding 2ir upon encircling the core of the defect. To understand 
the energies and the magnetic properties of such defects, it is useful to consider a London theory. Here we consider 
Phase 3 explicitly and state the results for the other phases. Allow {<j>i,4>2) to have phase windings of (n, m) times 
2tt respectively when encircling a defect. Taking (Aq x , Aq h , A_q x , A_q h ) = ^(e 1 ^ 1 , 0, e 1 ^ 2 , 0) /y/2, gives the London 
theory 

F=[ d 3 r^{^[(^! - 2eA t f + (c^ 2 - 2eA l f] + ±B?} (7) 

i 

where p StX — (k% + K2)\ip\ 2 , p s , y — («i — k^)!^ 2 , and p StZ — K3|V>| 2 - Closely related London theories appear in a 
variety of other contexts, including in a theory motivated by superconducting UPt3 [2(J, in two-gap superconductors 



2l| , and a description of fractional vortices (with non- Abelian vortex core states) in chiral superconductors [22| . Eq. [7] 
yields a supercurrent with components Ji oc p s ,i[di{4>\ + <fe) — 2e^4i]/2, far from the core of the defect, the minimum 
energy configuration has zero current and a contour integration then implies that the flux enclosed by a (n, m) defect 
is (n + m)$o/2. Consequently, the simplest defect encloses half a flux quantum. A similar analysis shows that Phase 
1 has only the usual superconducting vortices (the simplest containing flux $o) while Phases 2,3,4, and 5 all have 
defects that contain flux $q/2. Eq. [7] also shows that usual superconducting vortices have finite energy because the 
phase winding can be completely screened by the vector potential. The other defects have an energy that diverges as 
the logarithm of the system size. To help develop a physical picture of the fractional defects, it is useful to examine 
the induced CDW (and SDW) order near a <& /2 defect. Fig. Q] reveals that near such a defect in Phase 3 there is 
an edge dislocation in the CDW order when spatial uniformity of the defect along the z direction is imposed. If the 
defect is taken to be spatially uniform along the direction parallel to Q, then a screw dislocation appears in the CDW 
order. Note that the CDW order in this phase has half the periodicity of the PDW order, consequently, a dislocation 
in the CDW order can be interpreted as half a dislocation in the PDW order. The $0/2 defects of the other phases 
have a similar origin. This leads to the prediction that a PDW superconducting vortex containing a fraction of a flux 
quantum will be pinned to a dislocation in the accompanying CDW or SDW order. It is worthwhile pointing out that 
if Pi = in Eq.[TJ then Phases 4 and 5 will have <&o / 4 defects. A nonzero /J4 leads to a confinement potential between 
these $0/4 defects. However, it is still possible that the <&o / 4 defects play a role at short length scales, for example 
in determining the vortex core structure of integer flux superconducting vortices. 

We now turn to two physical consequences of the fractional flux defects. The first is in fluctuation driven vortex 
and dislocation physics in two dimensions. In this case, the fractional flux defects lead to non-superconducting 
CDW/SDW phases. The second consequence is that conventional superconducting vortices can decay into a bound 



FIG. 1: Fractional Vortex - Dislocation. Dislocation in CDW order accompanying a flux <E>o/2 defect in Phase 3. Since 
the periodicity of the PDW order is twice that of the CDW order, this corresponds to half a dislocation in the PDW order. 

state of fractional flux defects. This provides a mechanism for the appearance of CDW or SDW order inside a vortex 
core. 

The fact that single flux quanta vortices have short range interactions and fractional flux defects have long range 
interactions imply that novel vortex related physics can occur in two dimensions. This physics differs for the commen- 
surate and incommensurate cases and discussion of the commensurate case is left to the next paragraph. The theory 
for the incommensurate case follows from the London theory of Eq.[7]and resembles that for two band superconductors 
considered in Ref. (23|. To understand the behavior of the incommensurate case, consider initially taking the limit that 
the penetration depth A — * oo. In this limit the vector potential can be ignored, and a Bezerinskii Kosterlitz Thouless 
(BTK) transition, corresponding to the unbinding of $o/2 defects, occurs at a transition temperature T c * = |p s 
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24 , l25j |. In reality, since the penetration depth is finite, this will correspond to a crossover temperature at which 



the resistivity starts to fall. This is the case since superconducting vortices have a finite energy and therefore exist 
at any finite temperature. Note that this crossover temperature is half the value of the equivalent BKT transition 
temperature for a conventional superconductor with the same mean field T c and superfluid density. Consequently, the 
corresponding fluctuation regime will be larger. Now consider the case in which there exist thermally excited unbound 
superconducting vortex pairs. While such a state will not be superconducting, this does not imply that there is no 
longer any type of order present (as is the case for conventional superconductors). To see this, it is convenient to define 
9 = (</>i + </) 2 )/2 and <j) — 4>i ~ 4>2- The field 9 corresponds to conventional superconducting vortices and is disordered. 
The field (f> is uncoupled from the vector potential and therefore has vortex-like defects that are not screened. Vortices 
in this field correspond to dislocations in the CDW order. These dislocations cost a logarithmically divergent energy 
and therefore are not present at sufficiently low temperatures. This implies that there can be quasi-long range order 
in the CDW order parameter as defined through Eq. [4j The corresponding BKT transition at which the CDW order 
is removed occurs at T c d w < T* , above T c d w , strictly speaking, there is no order of any type. 

The commensurate case is more interesting since, in principle, it allows for the possibility to have T c d w > T* (where 
T* corresponds to a crossover temperature at which the resistivity starts to fall). To carry out a proper treatment, it 
is important to include the following term in the free energy density (which only exists in the commensurate case): 

e[A Q!c 4 (A_ Q /) 4 + A Q /(A_ Q /) 4 + A_ Qi 4 (A Qx *) 4 + A_ Qb 4 (A q /) 4 ]. (8) 

In the London theory for Phase 3, this leads to an interaction term cos(40i — 402). The renormalization g roup 
(RG) equations describing the CDW transition with such an interaction are the same as those derived in Ref. [26|], 
in which two-dimensional XY models subject to clock model-like symmetry breaking fields are considered. These 
RG equations imply that T c d w for the CDW order is enhanced relative to that of the incommensurate case. In the 
limit that Eq. [8] can be treated as a perturbation, we find that the enhancement of T c d w is given by T c d w {h,4) — 
T c dw(0) + ^15^2 l^ 4 l ^(^4) wnere ^4 — 2a 2 e|i/;| 8 and T c d w (0) corresponds to the BKT transition temperature with 
h 4 = Q. 

Finally, we turn to superconducting vortices containing a single flux quantum. These are the lowest energy vortices 
and they are created by magnetic fields. The key point is that a single $0 vortex can be either a conventional 
vortex or can be a bound state of fractional flux defects. This possibility is closely related to broken axisymmetric 
vortices discussed in the context of superfluid 3 He 2^, 28[, unconventional superconductors HH, and in FFLO 



superconductors [15| . The dissociation of a conventional vortex into a bound state of fractional flux defects provides 
a means to have a vortex core structure with CDW order. As a specific example, we have determined the structure 
of a vortex in a superconductor with non-vanishing d-wave and PDW superconductivity (the PDW order is that of 
Phase 3 or 4) . There exists a solution in which all components of the order parameter exhibit the same phase winding. 
Furthermore, the d-wave order parameter vanishes at the vortex core and the PDW order parameter is non-zero where 
the (i-wave order vanishes. This vortex will have CDW order (at the wavevectors listed in Table I) in the <i-wave 
vortex core as follows from Eq. |4j Such a solution may be relevant for understanding the observed CDW order inside 
the vortex cores of some underdoped cuprates [13] ■ 

In summary, we have examined the broken symmetry phases of PDW superconductors and have shown that the 
coexisting CDW or SDW order provides a means to distinguish between these phases. We have also shown that 
PDW superconductors exhibit topological defects that include fractional superconducting vortices which are coupled 
to dislocations in the coexisting CDW or SDW order. These defects can play an important role in stabilizing non- 
superconducting SDW or CDW phases. They are also important in understanding the physics of superconducting 
single flux quantum vortices where they can lead to CDW or SDW order inside vortex cores. Finally, in the context of 
the cuprates, we have examined the competition between PDW and translational invariant d-wave superconductivity 
and have shown that this prefers PDW phases with intrinsic CDW order and that an additional SDW or CDW appears 
due to the coexistence of d-wave and PDW superconductivity. 
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